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SPECTRAL ALBEDOS OF MiDLATITUDE SNOWPACKS 

Bhaskar Choudhury 

Hydrological Sciences Branch 
Goddard Space Fliglit Center 
Grcenbelt, Maryland 20771 (USA) 

ABSTRACT 

Spectral albedos of impure-ncnhomogeneous snowpacks, typical for midlatitudes, at wavelengths 
from 400 to 2200 nm are modeled through a numerical solution of the radiative transfer equation in 
the two-stream approximation. Discrete depth-dependent values of density, grain size and impurity 
concentration are used to characterize the s.nowpacks. The model is for diffuse incident radiation, 
and the numerical method is based on doubling and invariant imbedding. The effect of soot im- 
purities on snowpack albedos is illustrated: when a snowpack is semi-infinite (thickness of ten cen- 
timeters or more), soot reduces the albedos at visible wavelengths; however, for smaller snowpack 
thicknesses soot may increase the albedos at visible wavelengths. By adjusting soot contem and snow 
grain size, good quanti* f.ve agreement with some observations at the Cascade Mountains (Wash- 
ington) and at Point Barrow (Alaska) are obtained; however, the model grain sizes are found to 
be fifty to four hundred percent larger than the measured values. For satellite snowcover obser- 
vations, a model for effective albedo of partially snow-covered areas is developed and compared 
with some NOAA-2 observations of the southeastern United States. 
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SPECTRAL ALBEDOS OF MIDLATITUDE SNOWPACKS 


INTRODUCTION 

In mid-latitudes (30° - 60°), snowpacks accentuate their importance as a valuable water re- 
source for consumptive, industrial and agricultural use, and in flood control due to spring snow- 
melt runoff. Hydrologic importance of snowpacks, for example, in the Salt-Verde watershed in 
Arizona, and in the Sierra Nevada and the Cascade Range his been well documented (Rango, 

1975) The appearance and disippearance (except may be on high mountains) of snow at mid- 
latitudes creates a climatically sensitive transition zone because snow exerts considerable influence 
on weather and climate through numerous feedback processes (Hahn and Shukla, 1976; Kukla, 

1981 ; Reiter and Reiter, 1981 ). The snow data (e.g., areal coverage, depth, liquid water equivalent) 
from isolatcu weather stations are of limited use for hydrologic and climatic studies because they are 
point observations and extrapolation for areally representative values would be questionable. At- 
tempts are being made to obtain areally representative (including global) snow data, through remote 
sensing via visible, infrared and microwave radiometers (cf. McGinnis et al., 1975; Foster et al., 

1980; Matson and Wiesnct, 1981 ). 

Mid-latitude snowpacks are shallow (depth generally less than one meter except in high moun- 
tain environments), non-homogeneous (temperature gradient and melt-freeze metamorphisms and dep- 
osition of new snow on top of old metamorphosed snow) and impure (deposition of dust and soot 
from the atmosphere and scavenging of aerosols by the falling snow crystals). The objective of this 
paper is to formulate and study a model for spectral shortwave albedo of non-homogeneous shallow 
snoivpacks contaminated with impurities. Although incident radiation is assumed to be diffuse, the 
calculated albedos would be representative of those for the clear sky conditions at the solar elevation 
of about 35° (Wiscombe and Warren, 1980), i.e., about noon during the spring snow-melt period. 

In the next section the model is formulated in terms of a numerical solution of the radiative 
transfer equation. Illustrative results and comparisons with field and satellite observations are 
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then given. Some basic understanding of the radiative transfer process within a snowpack is gained 
through comparison with field observations, although interpretation of satellite data required address- 
ing additional issues. Previous studies (Warren and Wiscombe, 1 98J; Do/.ier et al., 1981; Choudhury 
et al., 1981) had concluded that absorbing impurities like soot would lower the visible albedo of a 
snowpack. The present calculations show that this effect of soot is not always true; when a snow- 
pack is not semi-infinite (i.e.. thickness less than 100mm, sec Wiscombe and Warren (1980)) the vis- 
ible albedo of an impure snowpack could actually be higher than that for pure snow. Comparison 
with satellite observations indicate that snew depth may be inferred from narrow band sensor data 
in the visible region, and this sensitivity to snow depth appears because the fractional snow covered 
area which determines the brightness of picture elements depends upon snow depth. 

FORMULATION FOR FIELD OBSFRVATIONS 

A snowpack is a delicate lattice whose nuclei are ice grains of various shape and size. The 
impurities could be attached to or imbedded in the ice grains or ligatures. In albedo studies 
(Bohren and Barkstrom, 1974; Choudhury and Chang, 1979; Wiscombe and Warren, 1980), a 
homogeneous snowpack of finite thickness has been modeled as isolated sphe' ical particles which 
scatter and absorb the ladiation incident on top of the snowpack. (An exception being Dunkle 
and Bevans (1956), who modeled the ice grains as facets.) The effect of absorbing impurities 
has been studied quantitatively within this schematic model of a snowpack (Warren and Wiscombe, 
1980; Choudhury et al., 1981>. \s an extension, a nonhomogeneous snowpack will be modeled as 
homogeneous layers (Fig. I), i.e., the model nonhomogeneous snowpack consists of layers which 
iiave characteristic values of density, grain size and impurity content. 

Within the above schematic representation of a nonhomogeneous snowpack, the albedo will 
be calculated through the two-stream approximation (Choudhury and Chang, 1979; Choudhury 
et al., 1981) of the radiative transfer equation. In this approximation, the radiation within the 
snowpack is divided into upwelling (F_) and downwelling (F^.) fluxes, which satisfy the equations 


dx 


(la) 


U) OJ 

* ( -F^. + y (I + g) F+ + y (I •• g) F.l 

dF. w cj 

“ — = I -F- + y (I +g) F_ + -y (I -g) F^) (lb) 

where the radiative transfer parameters y^, w, and g are called the extinction coefficient, the 
single scattering albedo and the asymmetry parameter of the scattering phase function. Tliese 
radiative transfer parameters vary with the wavelength and are functions of depth, x, througli 
their dependence on the physical properties of the snowpacic, namely, grain size, density and im- 
purity concentration. Calculation of the radiative transfer parameters from snowpack physical 
properties will be discussed before pa'senting the method of solution for eqns.(l). 

Since ice grains of a snowpack are two or more orders of magnitude larger than the wave- 
length, scattering by individual ice grains is strongly in the forward direction (the well known 
diffraction compoiu’U). In considering the radiative transfer process in powdered media, Hapke 
(1981) scrutinized the definition of the diffraction component as to whether this component 
should be treated as a part of the scattered radiation. As the inter-particle separation decreases, 
the particles enter into the geometric shadow of neighbouring particles, and Hapke (1981) con- 
cludes that when shadowing occurs the diffracted light must be treated as being indistinguishable 
from unscattered incident radiation. The condition for neglecting diffraction component as being 
a part of the scattered radiation is, for ice particles, given by (Hapke, 1981) 

< 0.94 (y)‘^^ (2) 

where Pj is density of ice, is snow density, r is the grain radius and X is the wavelength. The 
condition (2) implies that for natural snowpacks the diffracted radiation should almost always be 
neglected in calculating the radiative transfer parameters in eqns.(l). 

Noting that the extinction cross-section for particles two or more orders of magnitude larger 
than the wavelength can be taken as rrr- when the diffraction component is excluded, the extinction 


3 


cocfl’icicnl is given by (cl'.. Choudhury and Chang, I ‘>79) 


where pj is the density of ice (917 kg in’^), 


L 

4r 


(^) 

Pi 


( 3 ) 


The single scattering albedo is the ratio of scattering >nd extinction cioss-sections. This al- 
bedo inclusive and exclusive of the diffraction component, and co respectively, are related as 
(Van de Hulst, 1962) 


n = '/2 (1 + u;) (4) 

!n previous studies by the author and his co-workers (cf. Choudliury et al., 1981) the following 
equation obtained from Irvine and Pollack (1968) was used, 

CO exp (- 1.67 kr) (5) 

where k is the spectral absorption coefficient of the sphere. This ecpiation for a; lias the short- 
coming that for very large spheres (r -♦ <»), co tends to zero instead of approaching the correct 
value CO 2 i 0.066 (Irvine, 1965). This shortcoming arises because ecin.(5) neglects the rays re- 
Hected by the surface (Irvine and Pollack, 1968). Altiiough the error in calculating witli this 
shortcoming was found to be generally less than 5%, following Thomas (1952). the equation 
for CO chosen for this study is 

CO = 0.066 + 0.934 exp (- 1.75 kr) (6) 

A comparison of n obtained from (4) and (6) with the Mie theory results of Wiscornbe and 
Warren (1980) is given in Table I. 

The asymmetry parameter, G, inclusive of the diffraction component is related to the cor- 
responding parameter, g. excluding the diffraction cornpenent as (Van de Hulst, 1962) 
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1 + w 


( 7 ) 


Using the following parameterized equation for G (Choudhury et al., 1981), 

G « 0.886 u> + 0.978 (I - w) (8) 

one obtains 

0.022 

g - 0.886 - 0.092 w (O) 

<0 

A comparison of eqn. (8) with the Mie theory results of Wiscombe and Warren ( 1 980) is given in 
Table 2. Computationally, eqns. (6) and (9) are at least an order of magnitude less time consuming 
than the Mie theory. 

Although previous studies (e.g., Choudhury and Chang, 1979; Wiscombe and Warren, 1980; 
Warren and Wiscombe, 1980; Choudihury et al., 1981) on snow albedo had included the diffrac- 
tion component in calculating the radiative transfer parameters, it is shown in the Appendix that 
the calculated albedos from Choudhury and Chang’s (1979) model remain unchanged when the 
diffraction component is excluded, as it is done here. 

For pure snow, the absorption coefficients of ice from Grenfell and Perovich (1981) (from 
400 to 1400 nm) and Hobbs (1974) are used in eqn.(6). For soot contaminated snow, an effec- 
tive absorption coefficient for soot-ice mixture is calculated by averaging the absorption coeffic- 
ients of ice and soot by the corresponding volume fractions (see eqn.4 in Choudhury et al., 1981). 
The characteristics assigned to the soot particles are density 1 130 kg m'^, imaginary part of the 
refractive index 0.5 and particle volume 3.93 X 10“-* m^. The density and the refractive index 
of soot are for the assumed porosity of 50%(Chylek et al., 1981). (The results of Choudhury 
et al. (1981) were based on a higher value of soot density, 2260 kg m"^, which actually cor- 
responds to pure carbon. The porosity of soot is generally between 30 to 50%.) 
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For homoiBiMOUt tnowpackt (i.e.. p,. r and are independent of depth x) the lolution of 
eqni.(l) can be obtained analytically (Choudhury and Chang, 1979). A numerical method, how* 
ever, Its needed for nonhomogeneoui inowpacka. The numerical method uied in thii paper it 
baaed on doubling and invariant imbedding (Twoiney et al., 1966 and Grant and Hunt, 1969). 
Inatead of dlacuialng the numerical method in general terms, the basic steps in the algorithm im* 
plementation will be outlined. 

Step 1: initialization 

Hvlde anowpadc into NL homogeneous layers (Fig. ! ) and prescribe the I'^yer characteristics, 


Thickness : Ax j , Ax2 , Axj Axj^|^ 

Grain Radius : r|,r 2 ,r 3, 

Density : pj,P 2 ,P 3 Pnl 

Soot Content ; S| , 83 , S3 


After initialization execute the following steps for each wavelength. 

Step 2 ‘ Parameter Calculation 

Calculate the effective absorption coefficient (as discussed previously) and 


Single Scattering Albedo : c0| , 0)2. <03, , 

Extinction Coefficient : Tj . 72 » T3 7nl 

Asymmetry Parameter : g, , g2» g3 gjML 


Step 3: Doubling 

For each layer J (J ■ 1, 2, . . . NL), find an integer F such that 

Axi 

7j(— )<i0“^ (10) 

2 ^ 
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Cirfcuitte T|^(l) wmI Rn (I) throu|h the foUowliii recumnce relitioni 


TkO) 


TV|»> 

I-RV|W) 


(k - 1,2, 


. N) 
(II) 


Rfcd) - Rfc.i(J) ♦ 


t 2^.|(J) Rfe., (J) 


(12) 


with the lUrting relations 

Tod) ■ I “V^tj 


AXj 

II (I ♦g,)l — 





7j wj 0 - gj. 


2N 


(13) 

(14) 


Step 4: Invariant Imbedding 

Using NL values of T|,jd) and I^(J) (J ■ 1, 2, . . NL), execute the recurrence relations 

(L^^ 1.2 1) 


^ Tn(L+I)T(L) 

T(L+I) - s 

1-R(L)Rn(L+1) 

^ ^ t2n(L+1)S(L) 

RfL +1) • R(L) + 

I -R(L)R^(L+ 1) 


(15) 

( 16 ) 


with the starting relations 


T(I) - Tn(I) 

(17) 

R(I) - Rn (I) 

(18) 


Step y Spectral Albedo of Snow-Soil System 

If a is the spectral albedo of underlying soil, then the spectral albedo (A) of the snow-soil 
system is 

^ t2(NL)« 

A » R(NL) + (19) 

1 - a R(NL) 
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Spectral albedo of wet clay toll (Lulla, 1980) uaed In the calculation li ihown in Pig. 2. 

When the mow la cniited, -'me needi to conikler tte reflection due to the dielectric diicon" 
tinulty at top and bottom emit iurfa<»i md attenuation within the cruit layer. If the cruit layer 
la MHimed to be homogeneoui (l.e., volo of any trapped Impurltiei, air bubbles, etc.) with smooth 
surfaces (i.e., neidecting corrugation due to fused ice grains) then one can use Fresnel equations 
to calculate the reftectivities (Jerlov, 1976; Lyzenga, 1977), and an effective tiunsmission coefTic- 
lent (cf., Hottel e*^ il., 1968) for the layer. An equation analogous to (19) is obtained for crusted 
mow 

+ ^ { ( 20 ) 

I 1 - A ) 

where A b defined in eqn. 19, R^ Is Fremel reflectivity at the alr*ice interface for diffuse radia- 
tion 0.066; Jerlov, 1976) and Is elective transmisdon coefficient of crust layer for diffuse 
radiation, given by 

T. • (|.R^)(I-Vt}. . - - ' ■■ — { (21) 

I I -(r Rj.)2 ) 

where R|, Is Fresnel reflectivity at the ice-air interface for diffuse radiation (% 0.475; Lyzenga, 
1977) and r is the transmission function of the crust layer due to absorption, 

r ® exp ( ->/T k d) (22) 

k being the absorpuon coefficient and d the thickness of the crust layer. 

If the physical ^characteristics of a snowpack are known, the theory outlined above can be used 
to calculate its spectral albedos. In remote sensing via air- and satellite-borne sensors, the required 
physical characteristics are not generally known, and in fact a thrust of remote sensing is to infer the 
physical characteristics (e.g., the liquid water equivalent) from the sensor data. Also, high altitude 
observations would generally give the reflectivity of an area which is not uniformly covered with 
snow. The above theory needs to be amended for reflectively heterogeneous scenes. 
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In the next Mctlon, the it»ve theory will be compered with field obwrvetioni (uniformly enow* 
covered irea). The undentandii^ of the radiative tranifer procew gained through thia compariion 
will be used together with tome andllary data and an extenilon of the model to reflectively hetero> 
geneoui loenet for oomparii^ with latellite obtervationt. 

ILLUSTRATIVE RESULTS AND COMPARISONS 

Spectral albedo of deep aoot*amtaminated inowcover was studied previously by Choudhury 
et al. (1981). The study showed that the effect of soot is to tower the spectral albedo in visible 
and In neaNnf^ared up to about 1 100 nm. Furthermore, for the same wei^it fraction of soot, 
the reduction in albedo for a coarse n«ined snow was shown to be higher comparer! to a fine 
grained snow. Whereas the calculated spectral albedos for pure snow were generally higher than 
observations, the study substantiated Warren and Wiscombe*s (1980) work in providing a qualita* 
tive explanation for observed albedos. Choudhury et al. (1981), however, found some difficulty 
in matching concurrent observations of spectral albedos and flux extinction coefficients which 
they ascribed to be due to either approximations in the mode! or inaccuracies in the observations. 
Since tho present model is essentfally a modification of Choudhury et al. (1981) model to incor- 
porate layered nonhomogenities, the present model would not materially alter the conclusion! of 
Choudhury et al. when snow is deep and homogeneous. Choudhury et al. showed a comparison 
with the observations of Grenfell (1982). This comparison will be shown again, and will be com- 
plemented by comparisons with the observations of Grenfell et al. (1981) because the soot den- 
sity used by Choudhury et al. (namely, 2260 kg m*^) was a factor of two larger than the one 
used in this study. (1 have stated previously that the density 2260 ka m~^ is for s> lid carbon, 
and Warren and Wiscombe (1980) had correctly noted that soot being a porous carbon should 
have a lower density.) Some differences with Choudhury et al. also arise because the expression 
for the single scattering albedo and the absorption coefficient of ice in the spectral range 400 to 
1400 nm are slightly different in these two studies (see discussion in FORMULATION). 
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Compariion between ohieived (Grenfell, 1982; Grenfell et al., 1981) and calculated albedos 
for deep inowco%en are shown in Fip. 3(A, B, C). and observed and model physical parameters 
are given in Table 3. Observations in Fig. 3A are for new snow with ambient air temperature of 
0*C, while for Fip. 3B and 3C snow is old and ambient air temperatures are, respectively, *1 l*C 
and -f 13*C. All three figures show calculated albedos for pure (dotted curve) and soot contami* 
nated (solid line) r.ow, and it is seen that without absorbing impurities the calculated altn^do in 
the spectral range 400 to about 1000 nm are several percent higher than observations. In the 
spectral range 1500 to 2200 nm, the calculated albedos ore generally 3 to 10 percent lower than 
observations, which could be due to uncertainties in the absorption coefficient of ice (cf.. Wis> 
combe and Warren, 1980). The physical parameters given in Table 3 show that as the observed 
grain size increases, so does the model grain size. However, th^ model grain sizes are I.S to 2.5 
times larger than the observed values. 

Model calculations assume spherical grains, whereas natural snctv grains are rarely spherical. 
Since model calculation for nonspherical particles is more difficult than for spherical particles, it 
would be useful to find a measurable physical characteristic (e.g., (volume)* (surface area)*^^) 
which could be directly related to the model grain size (cf., O'Brien and Koh, 1981). A further 
point worth noting in Table 3 is the amount of soot needed to match the observations in Fip. 3 A 
and 3C which are from the same area (Cascade Mountains, Washin.^ton); the new snow in Fig. 

3A contains, according to the model, about six times more soot than old snow in Fig. 3C. Soot 
contamination in new snow would most likely be due to scavenging of atmospheric aerosols by 
the falling snow crystals. With aging one would expect the concentration of soot to increase due 
to fallout. Table 3, however, shows just the opposite trend for soot, although a substantial amount 
of dust was observed fo« Fig. 3C. One way to reconcile this discrepancy would be to postulate that 
snow>melt water rinses the micron size soot particles but not the larger dust particles. For matching 
the observations in Fig. 3B, Choudhury et al. (1981) used 0.5 ppmw (parts p-r million by weight) of 
soot, and the reason for a lower amount of soot in Table 3 is primarily the lower soot dnnsity used in 
this study. 
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Comparison with observed albedos for a two-layer deep snowpack is shown in Fig. 4, and 
the observed and model physical parameters are given in Table 3. Snow in the surface 30 mm is 
less than two days old, and the model grain size is somewhat unrealistically larger by about a 
factor of four than the observed size. Note also that soot in the surface layer Increases the al- 
bedo up to about 1000 nm. Fig. S(A, B) furthc. illustrates the effect of soot in shallow (less 
than a few centimeters deep) snowpacks. That absorbing impurities would increase the albedo 
when a snowpack is only a few centimeters deep is understandable. Choudhury et al. (1981) 
showed that the flux extinction coefficient, which determines whether a snowpack is semi-infinite 
at any wavelength, increases with soot content. The albedo of a few centimeters deep pure snow 
is lower in the visible because the incident radiation penetrates the snow layer, and thus is able to 
‘see’ a comparatively lower reflecting underlying soil (or metamorphosed snow). Soot increases 
the extinction coefficient and thus making a snow layer appear semi-infinite for thicknesses smaller 
than it would be for pure snow. W ith a thin (a few centimeUrs) layer of pure snow on a black sur- 
face, an interesting experiment would be to confirm that the albedo in the visible region to in- 
crease first and then decre&f-^ as absorbing impurities are added to the snowpack. An illustration 
of this interesting soot dep(..iuence of albedos is given in Fig. 6. 

Above comparisons showed that spectral albedos for uniformly snow covered areas can be 
matched if the model grain size and soot content are tre.^ted us adjustable parameters, and that 
fairly new snow contains ‘substaniial’ amount of absorbing impurities. In many observations by 
aircraft and satellite borne sensors the observed areas are not uniformly snow covered. Noting 
the relevance of satellite observed snow data in hydrologic and climatic studies, a semt-empirical 
analysis of data from the Very High Resolution Radiometer (VHRR) on board the NOAA-2 
satellite will be given. 

A linear model for sp''ctral albedo of a reflectively heterogeneous scene can be written as 

Aeff = S Aj (23> 

j® 1 
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where n is the number of components whose fractional areas visible to the sensor are fj and al- 
bedos Ai. (Note £ fi « l.> For a two-component scene (snow and vegetation), the equation 

J Ja| } 

Agff “ + (I ” 0 Aygg (24) 


would be 


where f is the fractional snow covered area visible to the sensor. 


Sat-Jllite observations are the bidirectional reflectance of the earth-atmosphere system gen- 
erally in narrow spectral bands. The albedos calculated here consider fluxes rather than in- 
tensities; however, 0*3rien and Munis* (1975) study shows a strong resemblance between the 
bidirectional reilectance and the albedo. Even with such resemblance, a direct quantitative 
comparison between the present model and the satellite observations can not be made betause 
of spectral absorption and scattering by atmospheric molecules and aerosols. Dozier et al.. 
(i981) have attempted to include the atmospheric effects, but a simpler approach <3f calcula- 
ting the ‘relative iciiectance* was followed by McGinnis et al. (1975). for the VHF,R data. A 
quantity equivalent to the relative reflectance will be the effective albedo normaliztrf by its 
maximum value. This normalization would minimize the atmospheric effects. 

For the VHRR sensor, bandwidth 0.6 - 0.7 /am, A^^g ^ 0.07 (Tucker, 1979; Kondratyev et 
al., 1981). The snow albedo does vary with snow thickness but, if comparison in Fig. 3A and 
illustration in Fig. 5 A are any guide, the fairly new snow observed by McGinnis et al. contained 
enough impurities to make two or three centimeters deep snow behave optically semi-infinite. 
From Fig. 3A we get Agnow * fractional snow covered area is known then eqn.(24) 

can be used to compute the effective albedo. 

The fractional snow-covered area that would be visible to a sensor depends, apart from sensor 
angle, on sun altitude, snow depth and surface features (e.g., canopy type, architecture and spac- 
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ing). With low sun altitude, the incident rays may get largely intercepted by the canopies, leav- 
ing the snowcovered ground mostly in shadows. Clearly, tall coniferous trees are more effective 
in 'hiding* the snow from ^!ie sensor view than coniferous shrubs and deciduous trees. 

In rugged terrains and low density coniferous forrests, the pixel brightness may increase In 
late winter or early spring (before snow melt) due solelv to higher solar altitudes. With increasing 
snow depth, furrows of a farm, scattered boulders, grass and other small plants and shrubs get 
covered with snow; only the taller plants protrude. The fractional snowcovered area will, there- 
fore, . 'crease as the snow depth increases, approaching the limiting value of 1.0 only if all objects 
within the pixel are covered with snow. 

A formal theory for fractional snow-covered area should certainly consider statistical or 
deterministic information about surface land features (e.g., height distribution and density of 
vegetation) within the pixel area. Instead of developing and verifying such a theory in this paper, 
I have used aircraft albedo observations of Kung et al. (1964) to calculate f. A brief discussion 
of Kung et al. experiment and calculation of f follows. 


A pair of upward and downward looking Kipp and Zonen solarimeters on board a Cessna 
310 aircraft were used by Kung et al. (1964) to study the spectrally integrated albedo over North 
America. The ratio of albedos for snow-covered and non snow-covered farm, mixed farm and 
forest, and forest areas in Wisconsin with a ground resolution of 21 x 70 meters is plotted as a 
function of snow depth in Kung et al. (1964). To analyze these plott?d data, a linear model for 
albedo would be 


+( 1-0 

^ ^veg^ ^ ^eg ^ 


( 25 ) 


where the angular brackets arc used to denote spectrally averaged values. 
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To calculate the fractional snow-covered area (0> it is thou^t that mixed farm and forest land 
data of Rung et al. would be somewhat morphologically compatible with 32 x 32 km squares VHRR 
data of McGinnis et al. (1975) for southeastern United States. The pertinent data are shown in Fig. 
7, and expressed parametrically as 


< Agff > 


60 + 3.6 X 


<Aveg> 


(26) 


60 + X 

where x is the snow depth (in mm). From eqns. (25) and (26) one can calculate f if snow and 
vegetation albedos are known. Rung et al. give the ..wan < Ayjg > for woody farm lands in Wis- 
consin as 0.144. Snow albedo is highly variable, and the factors affecting this albedo are meta- 
morphism, impurity contamination and snow depth. The snow in Rung et al. observations was 
at least two days old and, if the comparison shown in Fig. 3 A is any guide, contains enough im- 
purities t^ stipulate that a snow layer only two or three centimeters deep would be effectively 
semi-inflnite. If in Fig. 7 the largest ratio, 5.4 is assumed to be for f=l fi.e., uniform snow cover) 
then the corresponding snow albedo <c 0.778. The data in Fig. 7 spans the entire snow sea- 
son, and a snow albedo 0,778 is not totally unrealistic mean value for this data. Taking 
< > / < Aygg > =5.4, the fractional snow-covered area is obtained as (from Eqns. 25 and 26) 

0.6 X 


f = 


60 + X 


(27) 


From eqns. (24) and (27), the equation for pixel brightness (albedo) is 

4.2 + 0.6 X 

Aeff = (28) 

60 + X 

The effective albedo normalized to its value for x = 650 mm (i.e., A^fli = 0.56) calculated 
from eqn. (28) is shown in Fig. 8 together with the NOAA-2 observations of McGinnis et al. 
(1975). The calculated trend is in qualitative agreement with the observations; for a given snow 
depth, the calculated brightness value being somewhat higher than the observed mean value. 
McGinnis et al., however, had noted that the snow depth values that were associated with the 
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brightncM values may have been overeitimated. The effect of aiiociating a smaller snow depth 
value with the brightne»ses would be to shift the data points to the left and, if done, the agree- 
ment with calculation will improve. 

Note that the dependence of pixel brightness on snow depth is solely due to the dependence 
of the fractional snow-covered area on snow depth. Eqn. (27) for the fractional snow-covered 
area, obtained using Rung et al. (1964) data, implicitly includes sun angle and vegetation morpho- 
logy dependencies that was prevalent in Wisconsin, and its application to the NOAA-2 data for 
southeastern United States has not been rigorously justified. Comparison shown in Fig. 7 should, 
therefore, be tentative. If satellite observations are to be used to estimate snow depth, ab initio 
modeling of the fractional snow-covered area seems warranted. 

SUMMARY 

Midlatitudc snowpacks are a valuable water resource, and being in a climatic transition 
zone exerts considerable influence on weather and climate through numerous feedback processes. 
These snowpacks are generally nonhomogeneous and contain impurities, and an understanding of 
their spectral albedos is needed in energy balance calculations and remote sensing via visible and 
near-infrared radiometers. 

By characterizing a snowpack with depth-dependent density, grain size and impurity concen- 
tration, spectral albedos from 400 to 2200 nm are modeled through a numerical solution of the 
radiative transfer equation. The model is for diffuse incident radiation, and the numerical meth- 
od is based on doubling and invariant imbedding. The radiative transfer parameters (extinction 
coeffleient, single scattering albedo and asymmetry parameter of the phase function), which de- 
scribe scattering and absorption by the snow grains, are expressed analytically in terms of model 
grain size, measurable density and impurity content of the snowpack. The diffraction com- 
ponent of the scattered radiation is excluded in calculating the radiative transfer parameters to 
account for smaller inter-particle (grains) separation in snowpacks as compared to that in clouds 
or in the atmosphere. 
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Model calcuittioni ihow that when a snowpack if leveral centimeten deep, foot impuritief 
reducef the albedof up to about 1000 nm: however, for a few centimeterf deep fnowpackf, foot 
impuritief may increafe the albedo af compared to that for pure anow. Good quantitative agree- 
ment with fome obaeived spectral albedos at the Cascade Mountains (Washington) and Point 
Barrow (Alaska) are 'ihown, but the model ^ain sizes are larger than the measured values. The 
comparison points towards a better understanding of the meaning of the model grain size, as it 
is related to measurable dimensions of snow grains. 

For interpreting satellite snow-cover observations, a model for effective albedo of partially 
snow-covered areas is developed. An equation for the fractional snow-covered area, dependent 
upon the snowpack thickness, is obtained by analyzing a set of aircraft albedo observations. The 
predictions of the effective albedo model are in qualitative agreement with the NOAA-2 satellite 
data for the southeastern United States. It is concluded that remote sensing of snow depth using 
satellite data should address the question of fractional snow-covered area through ab Mtio model- 
ing. 
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APPENDIX: INVARIANCE OF CHOUDHURY-CHANG (1979) MODEL WITH RESPECT TO 


THE DIFFRACTION COMPONENT 


The aim of this appendix li to ihow that spectral albedo: and asymptotic flux extinction 
coefficients of snow obtained by Choudhury and Chang (*979) for diffuse incident radiation re 
main invariant with respect of inclusion or exclusion of the diffraction component in the radia* 
tive transfer parameters. 


Using the notation of this paper (FORMULATION) for the radiative transfer parameters, the 
spectral albedo (A) and the asymptotic flux extinction coefficient (0) obtained by Choudhury 
and Chang are 

(a Aoe- I) Aoo + (Aoo-a)e*2^ 

A ■ (A-I) 

(a Aoo - 1) + Aoo (Aoo - a) 


N (3 (1 - a) (1 - nc)l‘^ (A-2) 


a is soil albedo, h is snowpack thickness, is the extinction cross<section inclusive of diffrac* 
tion (2irr^), N is the number density of snow grains and A^o is the spectrai albedo for infinitely 
thick snowpack (h -*■ «), given by 


Aoo 


2(1 - O)^ 

1 

(1 - + (1 - OG)’^ 


(A-3) 


The radiative transfer parameters S2 and G inclusive of the diffraction component are related to 
the corresponding parameters excluding the diffraction component o) and g by relations (Van de 
Hulst. 1962) 

1 + wg 

G = (A-4) 

1 + w 


n = % (1 + w) 


(A-5) 
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One cm easily verify thet eqn. (A-3) for remeini invariant when eqns. (A'4) am* (A-S) 
are substituted in this equation. The eqn. (A*2) can also be seen to remain invariant when (A*4) 
and (A«S) are used, remembering that the extinction cross^ction when the diffraction 
component is excluded is 14 S^|. 

Thus the spectral albedos and the asymptotic flux extinction coefficients calculated from 
Choudhury and Chmg’s model will be identical to those calculated excluding the diffraction com- 
ponent 
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Table I. CompariscHi of co-albedos (I - O) from the parameterized equationfP) with the IGe results ini 
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Table 3. Observed and model physicd parameters for Figuna 3 and 4. 
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FIGURE CAPTIONS 


Fig. 1: Schematic iUuitration of the albedo model. Ax to layer thickneu, p to dentity. r to 

grain radiui, and S to the amount of toot in parti per million by weight. The lub- 
scripts are layer indices. F to the spectral incident flux and A is the spectral albedo 
of the snow>soil system. 

Fig. 2: Spectial albedo of wet clay soils. 

Fig. 3: Comparison with observed albedo for deep impure snowpacks. The dotted lines are 

for pure and solid lines are for soot contaminated snow. The physical properties of 
snow are given in Table 3. 

Fig. 4; Comparison with observed albedos for a twodayer snowpack. The dotted line is for 
pure snow in the top layer, the sub*layer being soot contaminated. 

Fig. S: Illustration of the effect of soot contamination for 20 mm thick snowpacks. (a) fairly 

new snow, (b) fairly old snow. The curves are labeled with the amount of soot (0.0. 

0.05 ppmw). 

Fig. 6: Spectral albedos from 400 to 1000 nm of a shallow snowpack for several discrete 

values of soot concentration. Spectral curves are labeled with the soot content 
(ppmw). This figure illustrates that soot first increases the albedos (left) and then 
reduces it (right). 

Fig- 7: The aircraft albedo data (o) and the parametric equation (dashed line) used for inter- 

fmg the fractional snow covered area. See text for the meaning of symbols. 

Fig: 8: The relative albedos (i.e., the albedos normalized to 100 percent at the snow thick- 

ness of 650 mm) for the NOAA-2 data (o). The dashed line is from Eqn. 28. 


26 




^*1 ' ^1 ' *"1 ' *1 


^*2 » ^2 ' *'2 ' *2 


SNOW 


SOIL 777777777777777777 ' 7777777777777777777777777777777777777 ‘ 


Figure 1. Schematic illustration of the albedo model. Ax is layer thickness, p is den- 
sity, r is grain radius, and s is the amount of soot in parts per million by 
weight. The subscripts are layer indices. F is the spectral incident flux and 
A is the spectral albedo of the snow-soil system. 
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Figure 2. Spectral albedo of wet clay soils. (Lulla, 1980) 
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Figure 3a. Comparison with observed albedo for deep impure snowpacks. TIk dotted 
lines are for pure and solid lines are for soot contaminated snow. The phys- 
ical properties of snow are given in Table 3. 
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SPECTRAL ALBEDO 



Figure 3b. Cc mparison with observed albedo for deep impure snowpacks. The dotted 
lines are for pure and solid lines are for soot contaminated snow. The phys- 
ical properties of snow are given in Table 3. 
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Figure 3c. Comparison with obser<^ed albedo for deep impure snowpacks. The dotted 
lines are for pure and solid lines are for soot contaminated snow. The phys- 
ical properties of snow are given in Table 3. 


31 



SPECTRAL ALBEDO 



Figure 4. Comparison with observed albedos for a two-layer snowpack. The dotted 

line is for pure snow in the top layer, the sub-layer being soot contaminated. 
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Figure 5a. Illustration of the effect of scot contamination for 20 mm thick snowpacks, 
fairly new snow. The curves are labeled with the amount of soot (0.0, 0.05 
ppmw). 
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Figure 5b. Illustration of the effect of soot contamination for 20 mm thick snowpacks, 
fairly old snow. The curves are labeled with the amount of soot (0.0, 0.05 
ppmw). 
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Figure 6. Spectral albedos from 400 to 1000 nm of a shallow snowpack 
for several discrete values of soot concentration. Spectral curves 
are labeled with the soot content (ppmw). This figure illustrates 
that soot First increases the albedos (left) and then reduces it 
(right). 
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Figure 7. The aircraft albedo data (o) and the parametric equation (dashed line) 
used for inferring the fractional snow covered are See text for the 
meaning of symbols. 


36 



RELATIVE ALREOO « 



Figure 8. The relative albedos (i.e., the albedos normalized to 100 percent at the snow 
thickness of 650 mm) for the NOAA*?. data (o). The dashed line is from 
Eqn. 28. 
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